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1 Introduction 



The calculus of variations is a beautiful and useful field of mathematics that deals with problems 
of determining extrema (maxima or minima) of functional 39 ,40, 58 . It starts with the simplest 
problem of finding a function extremizing (minimizing or maximizing) an integral 





J(y)= J F(t,y(t),y'(t))dt 



subject to boundary conditions y(a) = y a and y(b) — yb- In the literature many generalizations 
of this problem were proposed, including problems with multiple integrals, functionals containing 
higher-order derivatives, and functionals depending on several functions [35,37,45 . Of our interest 
is an extension proposed by Riewe in 1996-1997, where fractional derivatives (real or complex 
order) are introduced in the Lagrangian |50ii51j . 

During the last decade, fractional problems have increasingly attracted the attention of many 
researchers. As mentioned in [5], Science Watch of Thomson Reuters identified the subject as an 
Emerging Research Front area. Fractional derivatives are nonlocal operators and are historically 
applied in the study of nonlocal or time dependent processes [46] . The first and well established 
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application of fractional calculus in Physics was in the framework of anomalous diffusion, which 
is related to features observed in many physical systems. Here we can mention the report 38 
demonstrating that fractional equations works as a complementary tool in the description of 
anomalous transport processes. Within the fractional approach it is possible to include exter- 
nal fields in a straightforward manner. As a consequence, in a short period of time the list of 
applications expanded. Applications include chaotic dynamics [60j . material sciences [33], me- 
chanics of fractal and complex media [121132] , quantum mechanics [25]l3"T] , physical kinetics [61] , 
long-range dissipation 54 j, long-range interaction 1531 I55j . just to mention a few. One of the 
most remarkable applications of fractional calculus appears, however, in the fractional variational 
calculus, in the context of classical mechanics. Riewe [50,51] shows that a Lagrangian involv- 
ing fractional time derivatives leads to an equation of motion with nonconservative forces such 
as friction. It is a remarkable result since frictional and nonconservative forces are beyond the 
usual macroscopic variational treatment and, consequently, beyond the most advanced methods 
of classical mechanics [3D]. Riewe generalizes the usual variational calculus, by considering La- 
grangians that dependent on fractional derivatives, in order to deal with nonconservative forces. 
Recently, several different approaches have been developed to generalize the least action princi- 
ple and the Euler-Lagrange equations to include fractional derivatives. Results include problems 
depending on Caputo fractional derivatives, Riemann-Liouville fractional derivatives and oth- 
ers EMMIWMIMSMMMBHZMM- 

A more general unifying perspective to the subject is, however, possible, by considering frac- 
tional operators depending on general kernels [1 , 28 , 44 . In this work we follow such an approach, 
developing a generalized fractional calculus of variations. We consider very general problems, 
where the classical integrals are substituted by generalized fractional integrals, and the Lagrangians 
depend not only on classical derivatives but also on generalized fractional operators. Problems of 
the type considered here, for particular kernels, are important in Physics [T5]. Here we obtain gen- 
eral necessary optimality conditions, for several types of variational problems, which are valid for 
rather arbitrary operators and kernels. By choosing particular operators and kernels, one obtains 
the recent results available in the literature of Mathematical Physics [8] fT5Ul8] l24] . 

The paper is organized as follows. In Section[2]we introduce the generalized fractional operators 
and prove some of its basic properties. Section [3J is dedicated to prove integration by parts 
formulas for the generalized fractional operators. Such formulas are then used in later sections 
to prove necessary optimality conditions (Theorems [14] and [23]) . In Sections [4] [5] and [6] we study 
three important classes of generalized variational problems: we obtain fractional Euler-Lagrange 
conditions for the fundamental (Section [4]) and generalized isoperimetric problems (Section [6]), 
as well as fractional natural boundary conditions for generalized free-boundary value problems 
(Section[S]) . Finally, two illustrative examples are discussed in detail in Section!?] while applications 
to Physics are given in Section [8] in Section 18.11 we obtain the damped harmonic oscillator in 
quantum mechanics, in Section 18.21 we show how results from FALVA Physics can be obtained. 
We end with Section [9] of conclusion, pointing out an important direction of future research. 

2 Preliminaries 

In this section we present definitions and properties of generalized fractional operators. As par- 
ticular cases, by choosing appropriate kernels, these operators are reduced to standard fractional 
integrals and fractional derivatives. Other nonstandard kernels can also be considered as partic- 
ular cases. For more on the subject of generalized fractional calculus and applications, we refer 
the reader to the book [28\ Throughout the text, a denotes a real number between zero and one. 
Following [5], we use round brackets for the arguments of functions, and square brackets for the 
arguments of operators. By definition, an operator receives and returns a function. 



Definition 1 (Generalized fractional integral). The operator Kp is given by 



Kf[f}{x)=Kf[t^f{t)]{x)=p / k a (x,t)f(t)dt + q / k a (t,x)f(t)dt 



where P — (a,x,b,p,q) is the parameter set (p-set for brevity), x € [a, b], p,q are real numbers, 
and k a (x,t) is a kernel which may depend on a. The operator Kp is referred as the operator K 
(K-op for simplicity) of order a and p-set P , while Kp[f\ is called the operation K (or K-opn) 
of f °f order a and p-set P. 



Note that if we define 



G(x,t) 



pk a (x,t) if t < x, 
qk a (t,x) if t > x, 



then the operator Kp can be written in the form 



K a P [/] (x) =K a P [t^ f(t)} (x) = / G(x, t)f{t)dt. 

J a 

This is a particular case of one of the oldest and most respectable class of operators, so called 
Fredholm operators [231147] . 



Theorem 2 (cf. Example 6 of [23]). Let a € (0,1) and P = (a,x,b,p,q). If k a is a square 
integrable function on the square A = [a,b] x [a, b], then Kp : L2 ([a, b]) — > L% ([a, 6]) is well 
defined, linear, and bounded operator. 

Theorem 3. Let k a be a difference kernel, i.e., let k a 6 L\ ([a, 6]) with k a (x,t) — k a (x — t). 
Then, Kp : L\ ([a, b]) — > L\ ([a, b\) is a well defined bounded and linear operator. 

Proof. Obviously, the operator is linear. Let a G (0,1), P = (a,t,b,p,q), and / € L\{[a,b\). 
Define 

F(rA.-l \pk a (t~r)\-\f(T)\ iir<t 

^ T ^>-\ \qk a (r-t)\-\f(T)\ ifr>« 

for all (t, t) € A = [a, b] x [a, b]. Since F is measurable on the square A, we have 



I If F(r,t)dt\dr= f 

Ja \J a J Ja 



l/MI / \pk a (t-r)\dt+ \qk a (T-t)\dt 



dr 



< / l/(r)||b|-M|||fc a ||dT 

Hbl-klHlM- 11/11- 

It follows from Fubini's theorem that F is integrable on the square A. Moreover, 



ira/iii = 



t rb 

pl k a (t-T)f(T)dT + q k a (T-t)f( T )dT 



dt 



'■ /' (\P\ f\k a {t-T)\-\f{T)\dT + \q\ I \k a (r-t)\-\f(r)\dr)dt 



f ( J F{T,t)dr\ dt 



<\\p\-\q\[\\k a \\-\\f\\. 
Hence, K% : L t ([a,b]) -> L x ([a,b]) and \\K«\\ < ||p|-|q|| ■ ||fc ||. 



□ 



Remark 4. The K-op reduces to the left and the right Riemann-Liouville fractional integrals 
from a suitably chosen kernel k a (x,t) and p-set P. Let k a (x,t) — k a (x — t) = yT~i( x ~~ t) a ~ l : 

• if P = (o, x, b, 1, 0) , then 



K a P [/] (x) 



I» 



(x-tr-\f(t)dt=: a F x [f](x) 



is the standard left Riemann-Liouville fractional integral of f of order a; 
if P = (a, x, b, 0, 1), then 



K a P [/] (x) 



1 



I» 



(t - x) a -\f(t)dt =: x % [/] (x) 



is the standard right Riemann-Liouville fractional integral of f of order a. 

Corollary 5. Operators a L°^ lX L^ : L\ ([a, b]) — > L\ ([a, b]) are well defined, linear and bounded. 

The generalized fractional derivatives Ap and Bp are defined in terms of the generalized 
fractional integral K-oj>. 

Definition 6 (Generalized Riemann-Liouville fractional derivative). Let P be a given parameter 
set and < a < 1. The operator Ap is defined by Ap = D o Kp~ a , where D denotes the standard 
derivative operator, and is referred as the operator A (A-op) of order a and p-set P, while Ap[f], 
for a function f such that K p ~ a [f] £ AC ([a, b)), is called the operation A (A-opn) of f of order 
a and p-set P. 

Definition 7 (Generalized Caputo fractional derivative). Let P be a given parameter set and 
a G (0,1). The operator Bp is defined by Bp = K p ~ a o D, where D denotes the standard 
derivative operator, and is referred as the operator B (B-op) of order a and p-set P, while Bp[f], 
for a function f G AC ([a, b]), is called the operation B (B-opn) of f of order a and p-set P. 

Remark 8. The standard Riemann-Liouville and Caputo fractional derivatives are easily obtained 
from the generalized operators Ap and Bp, respectively. Let k\- a (x,t) = fci_ Q (x — t) = wrj ■, ■' 



• if P = (a, x, b, 1, 0) , then 
A a P [/] (x) = 



1 



-D 



( 



C ->• / K - t)- a f{t)dt 



(x) =: a D a x [/] (*) 



T(l -a)' 
is the standard left Riemann-Liouville fractional derivative of f of order a, while 

X 

1 



B a P [/] (x) 



(x-t)- a D[f}(t)dt=:^L^[f}(x) 



T(l - a) 

a 

is the standard left Caputo fractional derivative of f of order a; 
if P = (a, x, b, 0, 1), then 



-A P [f](x) 



1 



T(l -a) 



D 



(x) =: X I% [/] (*) 



s^ j(t-zr a f(t)dt 

is the standard right Riemann-Liouville fractional derivative of f of order a, while 

b 



-B a P [/] (*) 



•1 



(t-x)- a D[f](t)dt=:^D^[f](x) 



T(l -a) 

X 

is the standard right Caputo fractional derivative of f of order a. 



3 On generalized fractional integration by parts 

We now prove integration by parts formulas for generalized fractional operators. 

Theorem 9 (Fractional integration by parts for the K-op). Let a £ (0, 1), P = (a,t, b,p, q), k a be 
a square-integrable function on A = [a, b] x [a, b], and f,g G Li ([a, b]). The generalized fractional 
integral Kp satisfies the integration by parts formula 



b b 

g(x)K a P [/] (x)dx = f f(x)K%. [g] (x)dx, 



(1) 



where P* =< a,t, b,q,p >. 
Proof. Define 



F(r,t):= 



\pk a (t,r)\-\g(t)\-\f(T)\ iir<t 
\qk a (T,t)\-\ g (t)\-\f( T )\ ifr>t 



for all (r, t) G A. Applying Holder's inequality we obtain 



rb / ,-b 



IUits)*U-1 i/wi / 



l/MI / \pk a {t,T)\-\g{t)\dt+ \qk a (T,t)\-\g(t)\dt 



< 



< 



|/(V)| I / \ P k a (t,r)\-\g(t)\dt+ f \qk a (r,t)\-\g(t)\dt 



dT 



dr 



\f(r)\ 



2 I rb 



\pk a (t,r)\ dt 



\g(t)fdt 



2 I rb 



+ / \qk a (T,t)\ 2 dt 



\g(t)fdt 



dr. 



By Fubini's theorem, functions fc Q , r (i) := k a (t,r) and fc Q]T (i) := k a {T,t) belong to Li ([a, b]) for 
almost all re [a, 6]. Therefore, 



l/WI 



2 / cfc 



|pfc a (i,r)| eft 



lfl(*)l d« 



2 rb 



+ / |gfc Q (r,t)|^ / \g{t)\Ut 



dr 



hh l/WI lbWII 2 + 



< 



< 



2 / /■& 



IHM / l/WI ^ 



Qfca.r 



\pfca,T\ 



dr 



qk aT 



dr 



lffll 2 - ll/IMIIPM^ + lkMW < 00. 



Hence, we can use again Fubini's theorem to change the order of integration: 
6 b t b b 

g(t)KZim)dt = P Jg(t)dtJf(T)k a (t,T)dT + qJg(t)dtJf(T)k a (T,t)dT 

a a at 

b b b t 

= p [ f(r)dT ( g(t)k a (t,T)dt + q { f{r)dr I ' g{t)k a { T ,t)dt 



f(r)K^[g](T)dT. 



n 

Theorem 10. Let < a < 1 and P = (a, x, b,p, q) . If k a (x, t) = k a (x — t), k a , f £ L\ ([a, b]), 
and g € C ([a, b]), then the operator Kf, satisfies the integration by parts formula (fTj). 

Proof. Define 

\pk a (x-t)\-\g(x)\-\f(t)\ \£t<x 
\qk a (t-x)\-\g{x)\-\f{t)\ ift>l 

for all (t, x) £ A = [a, b] x [a, b]. Since g is a continuous function on [a, b], it is bounded on [a, b], 
i.e., there exists a real number C > such that \g(x)\ < C for all x £ [a, b]. Therefore, 



F(t,x) 



b I i-b \ r b 

F(t, x)dt dx - 



a \ ^ a 



< 



<c 



■b ft \ 

\f{t)\\ I \pk a (x-t)\-\g(x)\dx+ \qk a {t-x)\-\g(x)\dx 

J a j 

■ b rb > 

|/(*)| I / \pk a (x-t)\-\g(x)\dx+ \qk a (t-x)\-\g(x)\dx 



dt 



dt 



b pb \ 

\f{t)\\ I \pk a {x~t)\dx+ \qk a {t-x)\dx 



dt 



= C\\p\-\q\\\\k a \\\\f\\<^. 
Hence, we can use Fubini's theorem to change the order of integration in iterated integrals. □ 

Theorem 11 (Generalized fractional integration by parts). Let a £ (0,1) and P = (a,t,b,p,q). 
If functions /, KjZ a \g\ £ AC([a,b]), and we are in conditions to use formula ([T]) (Theorem^ or 



Theorem \10\) . then 
b 



g{x)B^[f]{x)dx= f{x)K x - a [g]{x)\ 



f(x)A a P , [g] (x)dx, 



(2) 



where P* =< a,t,b,q,p >. 

Proof. From Definition [7] we know that Bf,[f](x) = K x f a [D[f]] (x). It follows that 

!H.r)D)\[fU),l.r= j ! , [X )K 1 p - a [D[fH.rUh: 
J a 

By relation ([TJ 

b rb 

g(x)B a P [f]{x)dx = / D[f](x)K p - a [g](x)dx 

J a 

and the standard integration by parts formula implies @" 



{x)B a P [f](x)dx = f{x)K p Z a [g]{x)\ a - / f(x)D [K p Z a [g]] {x)dx 



U 



Corollary 12 (cf. [29:). Let < a < 1. If f,Jl~ a [g] £ AC([a,b]), then 



9 {x) c a D<* x [f]{x)dx=f{x)Jl- a \g]{x] 



I x—b 
I x—a 



f{x) x D« [g] {x)dx. 



4 The generalized fundamental variational problem 

By diF we denote the partial derivative of a function F with respect to its zth argument. We 
consider the problem of finding a function y = t H> y(t), t £ [a, b], that gives an extremum 
(minimum or maximum) to the functional 

J(y) = K% [t^F (t,y(t),y'(t),B^ [y] (*),#£ [y] (*))] (6) (3) 

when subject to the boundary conditions 

y(a) = 2/a, 2/(6) = 2/6, (4) 

where a, /3, 7 £ (0,1), -Pi =< a, b, b, 1,0 > and Pj =< a,t,b,pj,qj >, j = 2,3. For simplicity of 
notation we introduce the operator {-}pp defined by 

{yfp 2 \ (*) = (t,y(t),v'(t),B% [r » y(r)} (*),!% [r ^ y(r)] (*)) . 

With the new notation one can write (|3]) simply as J(y) = Kp F {y}p P (b). The operator Kp 

has kernel k a (x,t), and operators B p and if 7 , have kernels hi-p(t : r) and h 7 (t,r), respectively. 
In the sequel we assume that: 

(HI) Lagrangian F £ C 1 ([a, b] x R 4 ;R); 

(H2) functions A p , \t h> k a {b,T)diF {y} P ?p 3 (r) 



/T 7 



t ^ k a {b,r)d 5 F {yYpCpir) 



P,i 



D 



t i->- d%F {y}p p (t)k a (b,t) and t M> k a (b,t)d2F {y}p P (t) are continuous on (a, &); 



p 2 ,p 3 



1/3,7 



1-/3 



1/3,7 



(H3) functions i ->• ^Fjy}^ (t)k a (b,t), K p Z p |r h-> k a {b,r)d A F {y} P p 2 \ (r)J € AC([a,6]); 

(H4) kernels k a (x, t), h\^p(t, r) and /i 7 (£, r) are such that we are in conditions to use Theorems^ 

[I0]and[ll] 
Definition 13. A function y £ C l ([a, b]; R) is said to be admissible for the fractional variational 
problem ©H®, if functions Bp 2 [y\ and Kj, [y] exist and are continuous on the interval [a,b], and 
y satisfies the given boundary conditions (J3J). 

Theorem 14. If ' y is a solution to problem ([3|)-(|4]), then y satisfies the generalized Euler-Lagrange 
equation 



1/3,7 



d 



1/3,7 



k a (b,t)d 2 F{y} p «p 3 (t) - - {d 3 F{yy P2 ]p 3 (t)k a (b,t) 

- A p , [t ^ k a {b, t)8 a F {y}%\ (r)] (t) + K p , [r -> k a (b, r)d 5 F {y}%\ (r)] (t) = (5) 

for all t £ (a, b). 

Proof. Suppose that y is an extremizer of J . Consider the value of J at a nearby function 
y = y + sr], where e £ R is a small parameter, and r\ £ C 1 ([a, &]; R) is an arbitrary function with 
continuous B-op and K-op. We require that 77(a) = 11(b) = 0. Let 



J(y) = J(e) = J^ [t^F (t, y(t),y'(t),Bp: 2 [y] (t),Kp 3 [y] (t))] (6) 



d 

(ft 



k a (b,t)F ( t, y (*) + £77(i), ^ (y(t)+sr)(t)),Bp 2 [y + erj\ (t),K p [y + en] (t) ) dt. 



A necessary condition for y to be an extremizer is given by 
d2F{y}^P 3 V + d 3 F{y}%] P3 D[r ] } 



dJ 

de 



= 0<3-K£ i 



E=0 



d,F {y}^ P3 < [„] + d 5 F {yf P 2 P3 K\ [r,] 



> I I d ^ F {y}k\ (Mt)+d 3 F{y} P 2 P3 (t)- V (t) 



(6)=0 



(6) 



+ d 4 F {y}* 1 " (i)< [77] (t) + d 5 F {y}^„ (i)iH [77] (t) k a (b, t)dt = 



p 2 ,p 3 



Using classical and generalized fractional integration by parts formulas (Theorems |H [TO] and ITT|) , 



d< 



Wy}^ (i)fc a (MM*)[ - f j t (d 3 F{y}%\ (t)k a (b,t)) V (t)dt, 



^F{y}jC^(t)ka(J>,t)BP[r,](t)dt 



K 



1-/3 



r^k a (b,T)d 4 F{y}%\(r) (t)n(t) - A* T ^ k a (b,r)d A F {tf}&V (r) (t) »/(t)dt 



and 



MM) Wj/fc,, {t)lC P3 [rj\ (t)dt = / K},, \r ^ k a (b,r)d 5 F {y}^ (r) (*) r?(t)dt, 



where P* =< a, t, b, Qj 7 Pj >, j = 2, 3. Because 77(a) = 77(6) = 0, © simplifies to 



fca(M)fcF{y}£k (t) - j t (d 3 F{y}% P3 (i)fc a (M) 



At 



t h> k a (b,r)d 4 F{y}%\ (r)l (t) + #£. [r h- UM)Wy}& 3 (r) (t) }r)(t)dt = 0. 



We obtain §5$ by application of the fundamental lemma of the calculus of variations (see, e.g., (2"2l 
Section 2.2]). □ 

The next corollary gives an extension of the main result of |19j . 

Corollary 15. If y is a solution to the problem of minimizing or maximizing 

J(y) = J? [t ^F (t,y(t),y'(t), c a D? [y] (t))] (6) (7) 

in the class y € C 1 ([a, b];M) subject to the boundary conditions 

y{a) = Va, y(b) = y b , (8) 



where a, [3 E (0,1), F £ C 1 ([a,b] x IR 3 ;R) and r H> (6-t)"- 1 ^ (r,y(r),y'(r), £-D£ [y] (r)) /ias 
continuous Riemann-Liouville fractional derivative tD^ , then 

d 2 F (t, y(t), y'(t), c a Dl [y] (tj) ■ (b t)"' 1 - ± {dsF (t, y(t),y'(t), ?£>£ [y] (t)) ■ (b - t)^ 1 } 

+ t D% [r^(b-rr^d 4 F(r,y(r),y'(T), C a D^{y]{r))] (t) = (9) 

for all t £ (a, b). 

Proof. Choose k a (x,t) = YTa)( x ~ O" -1 ) hi-p(t, t) = r t-^_g\ (t — t)~@ , and P 2 =< a,t>&, 1,0 >. 
Then the A"-op, the ^4-op and the -B-op reduce to the left fractional integral, the left Riemann- 
Liouville and the left Caputo fractional derivatives, respectively. Therefore, problem (J?])-© is a 
particular case of problem ((3j)— ((4j) and ([9]) follows from ([5]) with d$F — 0. □ 

The following result is the Caputo analogous to the main result of [7] done for the Riemann- 
Liouville fractional derivative. 

Corollary 16. Let /3,7 € (0, 1). If y is a solution to the problem 



F 



(t,y{t),y\t), c a D%]{t), a F t [y]{t)) dt — > extr 



y e C 1 ([o,6];M) 
y(a) = y a , y{b) = y b , 



ther 



d 2 F (t,y(t),y'(t),CD P t {y](t), a my}(t)) - -d 3 F (t,y(t),y'(t),^[y](t), a r t [y}(t) 

'T*diF( T MT),AT),Zl%M(T),anMTJ)](t) 



tir b 



■t*b 



T <-> 



d 5 F(r,y(T),y'(r)^D^y}(r), a F T [y}(T))] (t) = (10) 



holds for all t G [a, b] . 



Proof. The Euler-Lagrange equation (|10l) follows from $5§ by choosing p-sets P\ —< a, b, b, 1, >, 
P 2 = P3 =< a,t, b, 1, >, and kernels k a (x, t) = 1, hx_p(t,T) = ph-m (^ — T )~^i an ^ ^(t, T ) = 

Remark 17. In i/ie particular case when the Lagrangian F of Corollary 1 1 6\ does not depend on the 
fractional integral and the classical derivative, one obtains from (|10[) the Euler-Lagrange equation 



5 Generalized free-boundary variational problems 

Assume now that in problem ©"(HI) the boundary conditions (0| are substituted by 

y(a) is free and y(b) = yb- (11) 

Theorem 18. Ify is a solution to the problem of extremizing functional (|3|) with (111[) as boundary 
conditions, then y satisfies the Euler-Lagrange equation ([5]) . Moreover, the extra natural boundary 
condition 



-1-/3 



dsF {y} p P2 \ (a)k a (b, a) + K l ~ p r -> 8 A F {j,}™ (r)k a (b, r) (a) = 



(8.7 
'2,-F-3 



(12) 



holds. 



Proof. Under the boundary conditions (|11[) . we do not require 77 in the proof of Theorem 1141 to 
vanish at t = a. Therefore, following the proof of Theorem 1141 we obtain 



WltfaTa WMMMaJ+rta)!^ r ^ 9 4 F{y}^ P3 (r)fc Q (6,r)j (a) 

d 



P2,P 3 

b 



+ J m id 2 F{y}%\ (t)k a (b,t) - - (d 3 F{y}%\ (t)k a (b,t)) 

t ^ d A F{y}^ P3 (r)fc a (6,r)l (*)+*!. [t •-» ^{y}^ (r)MM (t) I dt = II 



4 P 



(13) 



for every admissible 77 £ C l {[a, b];M) with 77(6) = 0. In particular, condition (fT3| holds for those 77 
that fulfill 77(a) = 0. Hence, by the fundamental lemma of the calculus of variations, equation ([5]) 
is satisfied. Now, let us return to ([T5)l and let 77 again be arbitrary at point t = a. Inserting J5]), 
we obtain the natural boundary condition (1121) . □ 



Corollary 19. Let J be the functional given by 

J(y) = a l%\t^F(t,y(t),atft[y)(t)) (//). 



Let y be a minimizer of J satisfying the boundary condition y(b) — y^. Then, y satisfies the 
Euler-Lagrange equation 

(b t) a - 1 d 2 F (t,y(t), °D?[y](t)) + tDS [r^(b- rJ-^F (r,y(r), ^y(r))] (t) = (14) 

and the natural boundary condition 

of [r^(b- T-r-^F (r, i/(t), ?Bfr(r))] (a) = 0. (15) 

Proof. Let functional © be such that it does not depend on the classical (integer) derivative 
y'{t) and on the If -op. If P 2 =< a,t,b, 1,0 >, hi-p(t — r) = ffi^al ft — r ) _/3 > an d fca(£ — t) = 
YT^j{ x — t) 01 ^ 1 ! then the i?-op reduces to the left fractional Caputo derivative and we deduce (fT4| 
and (|15j) from ([5]) and (TT21) . respectively. □ 

Corollary 20. Let J be the functional given by 

b 
J(y) = | F (t, y(t), y'(t), < [y](t), J% [y](t)) eft. 

a 

Lf y is a minimizer to J satisfying the boundary condition y(b) — yt,, then y satisfies the Euler- 
Lagrange equation 



d2F{ y y P2 \(t)-f t d 3 F{ y r P2 \(t)-A^ d,F{y} P 2 P3 (t)+K P * d 5 F {y}^ (t) = (16) 



and the natural boundary condition 

d 3 F {yf P2 \ (a) + K p / \d A F {yf"] (a) = 0. 



(17) 



Proof. Choose, in the problem defined by ([3]) and (TTT]) . k a (x,t) = 1. Then, equations (fTB")) and 
(H7J) follow from © and (JT2J), respectively. □ 
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6 Generalized isoperimetric problems 



Let (eR. Among all functions y : [a, b] — > K satisfying boundary conditions 

y{ a ) = Va, y(b) = S/b, 
and an isoperimetric constraint of the form 

I(y)=K$ l [G{y}' 3 £^](b)=t, 

we look for the one that extremizes (i.e., minimizes or maximizes) a functional 

J(y)=K£ 1 \F{yf£ p ](b). 



(18) 



(19) 



(20) 



Operators Kp , B p and KJ, , as well as function F, are the same as in problem (O - ®- Moreover, 
we assume that functional 1)1 9p satisfies hypotheses (H1)-(H4). 

Definition 21. A function y : [a, b] — > M. is said to be admissible for problem (|18p - (|20p if functions 
B p [y] and Kp [y] exist and are continuous on [a, b] , and y satisfies the given boundary conditions 
(|18p and the given isoperimetric constraint (|19p . 



Definition 22. An admissible function y £ C l ([a, &],R) is said to oe an extremal /or X if it 
satisfies the Euler-Lagrange equation §5§ associated with functional in (|T9"p . i.e., 

k a (b,t)d 2 G{y} p -2p 3 (t) j f (d 3 G{yf P 2p 3 (t)k a (b,t)) 

^ fc Q (6,r)a 4 G{y}^„ (r)l (t) + KL \t h-> k a (b,r)d 5 G {yfp 2 \ (r)l (t) = 0, 



A P , 



where P* =< a, t, 6, ^j,^ >, j — 2, 3, and t £ (a, o). 

Theorem 23. //y is a solution to the isoperimetric problem (|18p - (|20p and is noi an extremal for 
X, then there exists a real constant A such that 



k a (b,t)d 2 H{y} p ^p 3 (t) ± (d 3 H{y} p '2 Ps (t)k a (b,t)) 



A P' 



t ^ k a (b, r)d 4 H {y} P 2p 3 (r)l (t) + KL \t ^ k a (b, r)d 5 H {y}^ (r) (t) = (21) 



for allt £ (a, 6), where H(t, y,u,v,w) = F(t, y,u,v, w)—XG(t, y,u,v, w) and P* —< a, t, &, qj,Pj >. 

Proof. Consider a two-parameter family of the form y = y + £\r\\ + £2^2, where for each i £ {1, 2} 
we have T)i(a) = t]i{b) — 0. First we show that we can select £2^2 such that y satisfies (|19p. 

Consider the quantity X(y) = JCJ^ G{y}p P (0). Looking to Z(y) as a function of £i,£2, we 

define I{e\,E2) = ^(y) — £■ Thus, 1(0,0) = 0. On the other hand, applying integration by parts 
formulas (Theorems [5J QUI an d ITT1) , we obtain that 



9/ 

9£2 



(0,0) 



Ap, 



m (t)(k a (b,t)d 2 G{y} P 2p 3 (t) - j t [d,G{y} P 2 P3 (t)k a (b,t) 



t ^ k a (b, r)d 4 G {y} P Zp 3 (r)l (t) + /%. [r ^ fc a (6, r)d 5 G {y}^ (r) (t) ) dt, 



where P* =< a,t,b,qj,Pj >, j — 1,2. We assume that y is not an extremal for X. Hence, the 
fundamental lemma of the calculus of variations implies that there exists a function 772 such that 



11 



dl 

i:)62 



(0,0) 



7^ 0. According to the implicit function theorem, there exists a function £2(0 defined in a 



neighborhood of such that I(e\, £2(^1)) = 0. Let J(e 1; e 2 ) = J{y)- Function J has an extremum 
at (0,0) subject to 1(0,0) = 0, and we have proved that V/(0,0) 7^ 0. The Lagrange multiplier 
rule asserts that there exists a real number A such that V(J(0, 0) — A/(0, 0)) = 0. Because 



dJ 



(0,0) 



t (b,t)d 2 F{y}%\ (t) - ± (d 3 F{yf£ P3 (t)k a (b,t) 



A 1 ; 



\P,i 



T^k a {b,T)d A F{yY"{T) {t)+Kl r^k a (b,T)d 5 F{y} p ^ P3 (r) (t) ] m (t)dt 



and 



dl 



(0,0) 



k a (b,t)d 2 G{y}%] P3 (t) - j t (d 3 G{y}%\ (i)MM) 



A% 



k a (b,T)d 4 G{y} p ': p (T) (t)+K p , T^k a (b,T)d 5 G{y}%\(r) (t)) m (t)dt 



(8.7 
P 2 ,P 3 



one has 

b , 



J I k a (b,t)d 2 F{y}%\ (t) - j t (d 3 F{yf P 2 P3 (t)k a (b,t) 



A p 



T^k a (b,r)d 4 F{y} p P 2 P3 (r) (t) + K p , r ^ k a (b,r)d 5 F {y}™^ (r) (*) 
d 



\(k a (b,t)d 2 G{yf P 2 P3 (t) - - (d 3 G{y} P ^ Pa (t)k a (b,t)) 

r ^ k a (b,T)d 4 G{y} p ? P3 (r)l (t) + K P , [r * k a (b,r)d 5 G {y}%\ (r)l (t) 1 \ m {t)dt = 0. 



A fi 



From the fundamental lemma of the calculus of variations (see, e.g., [HJ Section 2.2]) it follows 

d 



Ub,t)d 2 F{y} P 2 P3 (t) - - (d 3 F{y} p ^ P3 (t)k a (b,t) 



A p , 



f P2,P 3 

d 



T^k a {b,T)d A F{y} P 2 P AT) (t)+K~L r^k a (b,T)d 5 F{y} P '7 P (T) (t) 



\[k a (b,t)d 2 G{yf P2 \ (t) - - (d 3 G{y} P ^ P3 (t)k a (b,t)) 

t h> k a (b,r)d 4 G{y} P 2 P3 (r)l (t) + KL fr » k a (b,r)d 5 G {y}^ (r) (t) I = n. 






that is, 

MM)Wv}&k (t) - I (a 3 F{y}^ P3 (t)*„(M)) 

r h+ fe a (6, r)5 4 ff M&V tol (*) + «!• r ^ fc Q (6, r)9 5 ff { y }™ (r) (t) = 



A p 



with # = F- AG. 



□ 
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Corollary 24. Let y be a minimizer to the isoperimetric problem 

J(y) =JS[t^F (t,y(t),?D?M(t))] (b) — > min, (22) 

!(?/) = „£ [i m. G (t lV (t), ^Df [v](*))] (6) - £, (23) 

2/(a) = y a , y{b) = y b . (24) 

If y is not an extremal of I, then there exists a constant A such that y satisfies 

(b^tr- 1 d 2 H(t, y (t)^Dny}(t))+tD b [r^(b-Tr- 1 d 3 H[r,y(r)^D^y}(r))](t)^O (25) 

for all t € (a, b), where H(t,y, v) — F(t,y,v) — XG(t,y,v). 

Proof. Let k a (x,t) = p^y(a; - t)"" 1 , h^^t.r) = r^^gj ft - r)^, Pi =< o, 6, 6, 1,0 > and 
P2 =< a,t,b, 1,0 >. Then the K-op and the -B-op reduce to the left fractional integral and the 
left fractional Caputo derivative, respectively. Therefore, problem ((22|) ((24)) is a particular case 
of problem (fl8 ]H pO ]) . and &$ follows from (J2TJ with d 3 H = d 5 H = 0. D 

Corollary 25. £ei y be a minimizer to 

J(y) = J F (t,y(t),y'(t), B p P2 [y](t),K], a [y}(t)) dt — > min, 
l(y)= J G(t,y(t),y'(t),B^[y](t),lC P3 {y}(t))dt = ^ 

J a 

y{a) = y a , y(b) = yb- 

If y is not an extremal of I, then there exists a constant A such that y satisfies 
d 2 H{y}^ P3 (t)-^d 3 H{y} P '2 P3 {t)-A P , [d 4 H {yf P2 \} (fi+Ify [d 5 H {y} P ' 2 \] (t) = (26) 

for all t € [a, b), where H(t, y, u, v, w) — F(t, y, u, V, w) — AG(i, y, u, v, w). 

Proof. Let in problem (JT5J) ([2TJ1) Pi =< a, b, 6, 1, > and kernel k a (x, t) = 1. Then, the generalized 
fractional integral K p becomes the classical integral and (|26|) follows from (f2l"j) . D 



7 Illustrative examples 

We illustrate our results through two examples with different kernels: one of a fundamental prob- 
lem ©-Q (Example l2"5|). the other an isoperimetric problem (flg ]) -(|2"0" ]) (Example [271) . 

Example 26. Let a, (3 E (0,1), £ e K, Pi =< 0,1,1,1,0 >, and P 2 =< 0,i, 1,1,0 >. Consider 
the following problem: 



J(y) = K a Pl 



t » tK^ [y] (t) + ^jl-(K^[y]{t)) 



(1) 



y(o) = i, y(i) = x 



rys(l-r) — 

(1 



-dr. 



r 2^2 



with kernel hp such that hp(t,r) = hp(t — r) and hp(Q) = 1. -Here i/ie resolvent rp{t) is related 



to the kernel hp(t) by rp(t) = C 



s l-> 



s/j^s) 



1 



(t), hf}(s) = C [t t— >• fya(t)] (s), where C and 



C 1 are the direct and the inverse Laplace operators, respectively. We apply Theorem \14\ with 
Lagrangian F given by F(t, y, u, v, w) = tw + vl — w 2 . Because 



y(t) 



1 



(l + £ 2 ) 5 ./o 



»*(*• 



1 



-dr 



(1+T 2 )- 



2\" 
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is the solution to the Volterra integral equation of first kind (see, e.g., Equation 16, p. 114 °f i4^0 



*k\M-¥g, 



it satisfies our generalized Euler-Lagrange equation ([3]), i.e. 

( 

t h> k a (b,r) 



K'i 



-K%[y](r) 



IfK'i^)) 2 ) 



(*) = o. 



In particular, for the kernel hp{t — r) = cosh(/3(£ — t)), the boundary conditions are y(0) = 1 and 
2/(1) = 1 + /3 2 (1 - y/2), and the solution is y(t) = l -^ + /3 2 (l - VT+1 2 ) (cf. [41, p. 22]). 

In the next example we make use of the Mittag-Leffler function of two parameters: if a, f3 > 0, 
then the Mittag-Leffler function is defined by 

E a>p {z)=Y J 



fc=0 



T(ak + P) 



This function appears naturally in the solution of fractional differential equations, as a general- 
ization of the exponential function |27| . 

Example 27. Let a, (3 <E (0, 1), ^gl, and £ ^ {^i} ■ Consider the following problem: 

J(y) = oi? 



Av) = oi? 



(y'+^iy})' 



y'+ c *n%'. 



(I) — > min, 



(27) 



2/(0) = , y (l) = £ £ x _ w (-(1 - r) 1 ^) Vl 4e 16e2 dr, 

which is an example of (|18|) - (|20[) with p- sets P\ —< 0,1,1,1,0 > and P2 =< 0,4,1,1,0 > and 
kernels k a (x — t) — rT^(x — t) a ^ 1 and h\^p{t — r) = T(1 _ g \ (t— t) • Function H of Theorem \23\ 

is given by H(t,y,u,v,w) = \J\ + (u + v) 2 — X(u + v) 2 . One can easily check (see \27\ p. 324]) 
that 



y(t)- 

• is not an extremal for X; 



E^i-it-rm^^dr 



(28) 



. satisfies y'+<<D?[y] = ^^. 
Moreover, (gSJ) satisfies (HTJ) /or A = 2£, i.e., 

/ / 

(l-tr^[y'(t)+^D^y](t) 



£ 

dt 



V 



\\ 



ty 



K \Jl + (y>{t)+cif t [y}{t)) 

( 

^(l-Tr- l (y'(T)+ZLt[y}(T)) 



4£ 



\ 



4£ 



v ^/l+(j/'W+ C -D?M(r)) 2 ) 

for all t 6 (0, 1). We conclude that (128j) is an extremal for problem ()27j) . 



(*) = 
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8 Applications to Physics 

If the functional ([3]) does not depend on B-op and K-op, then Theorem Q3] gives the following 
result: if y is a solution to the problem of extremizing 

J(y) = f F(t, y{t) lV '{t)) k a (b, t)dt (29) 



subject to y(a) = y a and y(b) = y^, where a € (0, 1), then 

d 2 F (t, y(t), y'(t)) - j f d 3 F (t, y(t),y'(t)) = — ^ ■ ^fc Q (6, t)d 3 F (t, y(t), y'(t)) . (30) 
We recognize on the right hand side of (|30|) the generalized weak dissipative parameter 

8.1 Quantum mechanics of the damped harmonic oscillator 

As a first application, let us consider kernel k a (b,t) — e a ( 6- *) and the Lagrangian 

L (y,y) = 2 m y 2 ~ v (y)> 

where V(y) is the potential energy and m stands for mass. The Euler-Lagrange equation (|30|) 
gives the following second order ordinary differential equation: 

y(t)-ay(t) = --V'(y(t)). (31) 

771 

Equation (|3Tj) coincides with (14) of [24], obtained by modification of Hamilton's principle. 

8.2 Fractional Action-Like Variational Approach (FALVA) 



We now extend some of the recent results of [15H18] , where the fractional action- like variational 
approach (FALVA) was proposed to model dynamical systems. FALVA functionals are particular 
cases of (f2"9"]l . where the fractional time integral introduces only one parameter a. Let us consider 
the Caldirola-Kanai Lagrangian [ MilfllTB] 



2 2 



L{t, V ,y)=m{t)[ V --u: 2 ^-), (32) 



which describes a dynamical oscillatory system with exponentially increasing time dependent mass, 
where ui is the frequency and m{t) = moe~~ fb e' yt = 777oe 7 *, 777,0 = 777oe _7b . Using our generalized 
FALVA Euler-Lagrange equation (j30|) with kernel k a (b,t) to Lagrangian (|32|) . we obtain 

y(t) + (6(t) + 7) y(t) + uj 2 y(t) = 0. (33) 

We study two particular kernels. 

1 . If we choose kernel 

MM) = ^-^pV+ 1 -i''+ 1 K, (34) 

defined in [55] , then the Euler-Lagrange equation is 

d 2 F (t, y(t),y'(t)) - f t d 3 F (t, y(t), y'(t)) = ( (1 '^^" ^ d 3 F (t, y(t), y'(t)) . (35) 
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In particular, when p — > 0, (|34l) becomes the kernel of the Riemann-Liouville fractional 
integral, and equation (j3"S"j) gives 

d 2 F (i, y(t), y'(t)) - j t d,F (t, y(t),y'(t)) = ^-| d 3 F (t, y(t), y'(t)) , 

which is the Euler-Lagrange equation proved in [16J . For p =^ 0, we have 

v ; 6P+ 1 - tP+ l 

Therefore, both at the very early time and at very large time, dissipation disappears. More- 
over, if p — s- 0, then 

1 — a [ ift^oo 



S(t) = 



b-t [i^ ift->0. 



This shows that at the origin of time, the time-dependent dissipation becomes stationary, 
and that at very large time no dissipation, of any kind, exists. 

2. If we choose kernel k a (b,t) = (cosh 6 — coshi)" , then 

d 2 F (t, y(t), y'(t)) - -d 3 F (t, y(t), y'(t)) =-(a- l)—^—-^-d 3 F (t, y(t), y'(t)) (36) 



and 



... . . sinh t [a — 1 if £ — > oo 

6(t) = -{a - 1) — > < 

w v ; coshfe-cosh£ |0 if t -)• 0. 



In contrast with previous case, item 1, here dissipation does not disappear at late-time 
dynamics. 



We note that there is a small inconsistence in [16], regarding to the coefficient of y(t) in (|33|) . 
and a small inconsistence in |18) . regarding a sign of (|36|l . 

9 Conclusion 

In this article we unify, subsume and significantly extend the necessary optimality conditions 
available in the literature of the fractional calculus of variations. It should be mentioned, however, 
that since fractional operators are nonlocal, it can be extremely challenging to find analytical 
solutions to fractional problems of the calculus of variations and, in many cases, solutions may not 
exist. In our paper we give two examples with analytic solutions, and many more can be found 
borrowing different kernels from the book [47]. On the other hand, one can easily choose examples 
for which the fractional Euler-Lagrange differential equations are hard to solve, and in that case 
one needs to use numerical methods [2 , 6 . 48 , 49 . The question of existence of solutions to fractional 
variational problems is a complete open area of research. This needs attention. Indeed, in the 
absence of existence, the necessary conditions for extremality are vacuous: one cannot characterize 
an entity that does not exist in the first place. For solving a problem of the fractional calculus 
of variations one should proceed along the following three steps: (i) first, prove that a solution 
to the problem exists; (ii) second, verify the applicability of necessary optimality conditions; (iii) 
finally, apply the necessary conditions which identify the extremals (the candidates). Further 
elimination, if necessary, identifies the minimizer(s) of the problem. All three steps in the above 
procedure are crucial. As mentioned by Young in [59] , the calculus of variations has born from 
the study of necessary optimality conditions, but any such theory is "naive" until the existence of 
minimizers is verified. The process leading to the existence theorems was introduced by Leonida 
Tonelli in 1915 by the so-called direct method [56) . During two centuries, mathematicians were 
developing "the naive approach to the calculus of variations" . There was, of course, good reasons 
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why the existence problem was only solved in the beginning of XX century, two hundred years after 
necessary optimality conditions began to be studied: see [13U57J and references therein. Similar 
situation happens now with the fractional calculus of variations: the subject is only fifteen years 
old, and is still in the "naive period" . We believe time has come to address the existence question, 
and this will be considered in a forthcoming paper. 
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